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ABSTRACT 



Context. The study of the linear and circular polarization in AGN allows one to gain detailed information about the properties of the 
magnetic fields in these objects. However, especially the observation of circular polarization (CP) with single-dish radio-telescopes 
is usually difficult because of the weak signals to be expected. Normally CP is derived as the (small) difference of two large numbers 
(LHC and RHC); hence an accurate calibration is absolutely necessary. 

Aims. Our aim is to improve the calibration accuracy to include the Stokes parameter V in the common single-dish polarimetric 
measurements, allowing a full Stokes study of the source under examination. 

Methods. A detailed study, up to the 2nd order, of the Miiller matrix elements in terms of cross-talk components allows us to reach 
the accuracy necessary to study circular polarization. 

Results. The new calibration method has been applied to data taken at the 100-m Effelsberg radio-telescope during regular test 
observations of extragalactic sources at 2.8, 3.6, 6 and 11 cm. The D-terms in phase and amplitude appear very stable with time and 
the few known values of circular polarization have been confirmed. 

Conclusions. It is shown that, whenever a classical receiver and a multiplying polarimeter are available, the proposed calibration 
scheme allows one to include Stokes V in standard single-dish polarimetric observations as difference of two native circular outputs. 

Key words. Polarization - Radio continuum: galaxies - Techniques: polarimetric - Methods: observational 
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1. Introduction 

The study of circular polarization (CP) and accurate multi- 
frequency observations of linear polarization (LP) provide 
probes of AGN magnetic field properties. CP, LP and spec- 
tral information can be used to constr ain the low ene rgy end 
of the relativistic particle distribution IBecker ti l2003ll. derive 
magnetic field order, strength and geometry [Gabuzda, 2006] 
and make inference s about the particle content of the AGN jets 
llWardle et al.lll998Tl . So far several studies aimed at studying CP 
from AG N have been carried out with interferometers: ATCA at 
5 GHz iRaroeret aUlIOOOll . VLBA at 5.0 GHz a nd 15.0 GHz 
llHoman & WardleL 1 1999(1 and the VLA at 8.4 GHz llBower et all 
2002]. On the other hand for instance t he Michigan 26- metre 
radio-telescope at 4.8 GHz and 8.0 GHz llAller et all I2003H has 
been used to perform single-dish CP observations. 

Stokes V, which measures CP through the ratio V/I, has of- 
ten been discarded from the polarimetric observation strategy 
and calibration, due to its typically low values, leading to a 
strong preference for radio-astronomical receivers to provide cir- 
cularly polarized outputs, as this configuration allows efficient 
measurement of Stokes /, Q and U. 

When a single-dish telescope delivers two opposite circular 
outputs, it is relatively easy to retrieve information about Stokes 
/ (sum of the outputs) as well as Q and U (cross-correlation 
of the outputs), but the estimate of Stokes V (difference of the 
outputs) can be severely contaminated by the receiver internal 
noise and gain fluctuations and by small intrinsic asymmetries in 
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the receiver channels structure (e.g. differences between cables 
lengths, different internal gains, unwanted internal reflections). 
Some of these effects cause a small percentage of contamination 
between the signals belonging to the left and the right channel 
(D-terms), in amplitude and phase. 

To fill the gap that Stokes V leaves in the single-dish po- 
larimetric observations, the most common receiving architec- 
ture has been examined. A full Stokes calibration scheme has 
been derived by identifying at each step the possible contamina- 
tion that may be introduced by the devices used in the receiving 
chain, along with the dependance of the 4x4 Miiller matrix ele- 
ments on the D-terms, up to 2nd order. The result has been ap- 
plied to several observations carried out at the Effelsberg 100-m 
telescope. 

The application of an accurate 4x4 Miiller matrix to the 
observed raw data offers the possibility of performing multi- 
frequency full Stokes observations of the sources, giving a com- 
plete description of their polarization spectra. 



2. Basics of polarimetry 

2. 1 . Polarized waves 

By focusing on the electric field only (the magnetic field can al- 
ways be derived from E = -c x B) a wave that is travelling 
along the z-axis with complex amplitude Eq can always be de- 
composed into two components, one along the jt-axis and one 
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along the y-axis, as follows: 

~E = iE x e i(kz -" ,+s ' ) + }E y J braM > ) 

= [iE x e i6 ' + JEy^f] e i(kz -^ ( 1 ) 

= E e i(kz - M) 

where i and / are unit vectors in the x and y direction. 

Polarization is a phenomenon which can occur in transver- 
sal waves and is described by the "trace" left by a wave pass- 
ing through a plane perpendicular to the propagation direction. 
A well defined polarization state (linear, circular, elliptical, also 
called "mode") can be recognized in all those cases where a con- 
stant relation between E x , E y and 6 X , 6 y exists and reflects the 
existence of a defined ordered motion of the electric vector. 

Any polarized wave can be considered as composed of two 
linearly polarized waves, as shown in Eq. ([TJ, but any pair of or- 
thogonally polarized waves (e.g. left circular/right circular, left 
elliptical/right elliptical with axis at 90°, etc.) can be used as a 
basis for description of a polarization state. The most convenient 
choice is typically the polarizations to which the telescope out- 
puts respond. 

2.2. Jones vector and Muller matrix in the LR circular 
representation 

The polarization properties of the radiation can be described us- 
ing different vectors. The Jones vector describes the information 
related to amplitude and phase of the observed signal compo- 
nents, thus it is the preferred one to describe mathematically the 
receiver internal devices. The Stokes vector reflects the relation 
between the signal components and defines the overall polariza- 
tion of the signal, thus it is the preferred one to describe the 
source physical characteristics. 

The Jones and Stokes vectors carry the same information 
and are linked by a set of transformation equations. As a rule 
of thumb, the basis for the transformation should be chosen ac- 
cording to the polarization of the outputs delivered by the tele- 
scope under examination. The following equations are related to 
the circular L R fr ame, which w e adopted in this paper following 
iKrauslll 198611 and lHeilesI |l2002tl . 

The Jones vector, J , and the Stokes vector, S , are defined as 
follows 



J = 



E L 




E L e- iSt 


Er 




E R e- i5, < 



i = {e 2 l ) + (e 2 r ) = eie l + e r e r 

Q = 2 (E L E R cos A> = E L E R + E R E* L 
[7 = 2 (E L E R sin A) = i (e r E* l - E L E R ) 
V = {e 2 l )-{e 2 r )=E* l E l -E r E r 



(2) 

(3) 

(4a) 
(4b) 
(4c) 
(4d) 



The Muller matrix relates the polarimetric properties of the 
source, i.e. the Stokes parameters, of the source with the mea- 
surements from the polarimeter at the telescope. Hence it de- 
scribes the instrument polarization characteristics. An ideal in- 
strument that measures Stokes parameters perfectly, would be 

— > — > 

characterized by a unitary Muller matrix S m = MS , 



M 



'run W13 mu 

m 2 i m 22 m 2 3 m 2 4 

OT31 OT32 »?33 »t34 

WI41 W42 7W43 m44 J 



(5) 



(6a) 
(6b) 
(6c) 
(6d) 



I m = m\il s + mizQs + mi 3 U s + m u V s 
Qm = mzih + m 22 Q s + m 23 U s + m 2A V s 
U m = m 3 il s + m-szQs + m 33 [7 s + m iA V s 
V m — m A \I s + m A2 Q s + nii^U s + m AA V s 

where I m , Q m , U,„ and V m are the measured Stokes parameters 
and 7 S , Q s , U s and V s are the true Stokes parameters describing 
the source polarization. 

When handling linear polarization, the polarization angle is 
measured with respect to some instrumental zero point, usually 
related to a reference direction on the celestial sphere. During the 
source tracking, the reference system of an alt-azimuth mounted 
telescope rotates with respect to the source, that is, the parallac- 
tic angle of the sources (and consequently i ts polarization an- 
gle) rotates during the observation. Following lTurlo et al.l [l985] 
the time -dependent rotation component B of the Muller matrix 
can be extracted and applied separately from the instrumental 
component of the Muller matrix, T is assumed to be constant 
with time. The overall Muller matrix can be then expressed as 
M = TB, where B is the rotation matrix of a linear rotating sys- 
tem (rotation of ■&) 



B = 



1 0\ 
cos 2§ - sin 2i? 
sin2i? cos2# 

00 01; 



B = 



1 ^ 
cos 2$ sin 2§ 
- sin 2d- cos 2§ 
1 

(7) 



where the asterisk stands for complex conjugate and A = 5 r -5l- 



3. Application to a single-dish observation 

iMcKinno n 1 1992] and iHamaker et al.l Jl996ll published major 
works o n how to manag e polarimetry with a radio interferom- 
eter and iJohnst on 1 2002] adapted their study to single-dish an- 
tennas equipped with linear dipoles. Nowadays telescopes are 
commonly equipped with scalar feeds and hybrid or wave-guide 
elements that supply circularly polarized outputs. In the follow- 
ing we describe, step by step, how this instrumentation can af- 
fect, in terms of instrumental polarization, the measurement of 
the incoming radiation. A typical receiving chain is sketched in 
Fig- HI In Table[T]we summarize our notation. 

3.1. The antenna 

The targets of single-dish observations are typically point-like 
with respect to the antenna beam and the observables are the 
integrated Stokes parameters and the full Stokes spectra of the 
source. The local CP and LP can be observed only resolving the 
source with interferometric techniques, and are often averaged 
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Fig. 1. Receiver sketch. The block diagram shows the main parts of a typical radio-astronomical receiver: the scalar feed, sensitive 
to all the polarization states; the directional coupler, that couples the linearly polarized noise diode signal into the source radiation; 
the polarizer, that transforms the signal into the circularly polarized representations; the OMT, that splits the orthogonally polarized 
signals (polarizer and OMT are sometime realized as a single device); the low noise amplifiers (LNA); the down-conversion parts, 
here represented in a single block (labeled "conversions"); the total power detectors (TP); the two multipliers, here shown as a single 
device, that take as input the LHC and RHC channels to give out Stokes Q, and the LHC and the 90° phase shifted RHC to give out 
U. The picture shows also equations giving the electric field strength, or voltage in the parts where propagation is as current on a 
wire, how those depend on the gain and the D-terms factors and how they are combined during each propagation step, see Ch. 4. 1 
and Ch. 4.2 for a complete derivation. 
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Table 1. Here the notation adopted to refer to different quantities 
is summarized. All complex quantities are expressed in the form 
A = Ae' p , where A is the amplitude and p is the phase. P is the 
phase difference between the two channels. 
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out when the source is not spatially resolved (beam depolariza- 
tion). 

An ideal rotationally symmetric telescope does not change 
the state of polarization of the incident radiation (though circu- 
lar polarization changes its state every ti me it is reflected, due 
to cha nging direction of propagation) see McGuire & Chi pman 
lfl99l . 

The first deviation from the ideal case is due to the blockage 
by the subreflector and its supporting structure. The radiation 
scattering causes a cross-polarization effect, both in amplitude 



and phase, which can be widely reduced by the usage of an ap- 
propriate structure design. Coupling the antenna with corrugated 
horns or dual mode horns with circularly symmetric patterns of- 
fers cross-polarization efficiencies closer to 1 . 

The parabolic shape is usually optimized at a well defined 
elevation, typically around 45°. When the antenna observes far 
from this value, its shape deviates from the original parabola 
(and the symmetry could be broken), unless it is equipped with 
some correction mechanism that compensates the gravitational 
deformations (active surfaces or homology). As an example of 
the deformation involved, Table [2] lists the deviation from the 
best-fit parabola for a typical dish without compensating tech- 
niques. 



Table 2. Deviation from the best fit parabolic surface, typical 
values 



Elevation 


Deviation 


Deviation 


[°] 


[rms mm] 


[% @ 6 cm] 


30 to 60 


0.5 


0.8 


20 to 75 


0.8 


1.3 
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3.2. The feed 

The feed is used to obtain a Gaussian beam pattern and con- 
verts the electromagnetic field from propagating in free space to 
propagating in waveguide. The most common type, coupled to 
Cassegrain and Gregorian antennas, is a conical corrugated horn 
(often called "scalar feed"), as it provides a Gaussian beam pro- 
file and is characterized by negligible cross-polarization. Unlike 
linear dipoles or helical antennas, the scalar feed usually does 
not have a single defined polarization mode, as it can support 
linear, circular and dual polarization, depending on its design 
and on the electronic devices that are coupled to it. A typical 
scalar feed provides a single output: some additional devices are 
needed to split it in the polarized components. 

A lack of symmetry in the feed conical shape could lead to 
some degree of polarization aberrations. As an example of devi- 
ation from the ideal feed shape, the typical mechanical accuracy 
is 0.02 mm to 0.03 mm. This would lead, at 6 cm (nearly 600 mm 
of diameter, in secondary focus) to a possible non-circularity of 
0.003% to 0.005%. 

For off-axis feeds (multibeam or focal plane receivers that 
permanently hosts several different feeds) the rotational sym- 
metry is broken. As a result, a certain amount of instrumental 
polarization enters and increases the polarized signal I p and the 
antenna patterns for the two orthogonal polarization channels of 
the receiver are different (beam squint). The result is a constant 
pointing offset between the two opposite beams w hich depends 
on th e distance of the feed from the optical axis llFiebig et al.L 

mm. 

3.3. The directional coupler for calibration 

A radio-astronomical receiver suffers from a certain amount of 
gain instability, due to electronic noise added by each compo- 
nent, thus two measurements of an identical power source could 
yield slightly different voltage values. To stabilize this fluctua- 
tion an internal noise of known temperature (Teal) is added to 
the astronomical signal using a directional coupler. This calibra- 
tion signal is linearly polarized and is injected into the receiver 
with a known polarization angle. The Teal is also used to convert 
the measured voltage into kelvin. 

The value of the Teal is usually known with an accuracy of 
3% to 5% over a wide bandwidth (e.g. 500 MHz for the 6 cm 
receiver at the Effelsberg telescope). 

A typical value of the stability is of 0.01 dB/K, this number 
must be interpreted as a Teal variation of 0.2% per each 1 K 
variation in the receiver room temperature (which is usually kept 
constant). Assuming a Teal nominal value of 2 K and a thermal 
excursion inside the room of + 5 K a value of (2.00+0.02)K is 
obtained. 

3.4. The polarizer and the ortho-mode transducer (OMT) 

The polarizer and the OMT split the signal in two orthogonal cir- 
cular polarized components. This process is not perfect and some 
fraction of one component can affect the other component. This 
effect is known as "cross-talk" and is quantified by the so called 
"D-terms" (values < 1% are usually considered very good). 

The ensemble of feed, polarizer and OMT is completely 
analogous to a pair of circularly polarized feeds, although the 
instrumental effects on the polarization measurements are differ- 
ent. 

At this point in the receiving chain the following can be ex- 
pected: 



• Antenna: the non-perfect rotational symmetry can cause 
diattenuation and birefringent effects (that is, some instru- 
mental polarization and some conversion among Q, U, V). 
For circular polarization, a change in sign must be expected 
at each reflection. These effects are usually cancelled by 
appropriate technical equipment (active surface, appropriate 
cabling to compensate the change in sign) and by a calibration 
that compensates possible residual components (e.g. antenna 
gain curve as function of elevation). 

• Feed: small mechanical imperfections cause additional 
diattenuation and birefringent effects, but at low frequencies 
this should be a negligible contribution, as, e.g. at 6 cm, the 
imperfections are typically of the order of 0.01%. 

• Directional coupler: the phase of the transmitted signal is 
not constant over the bandwidth, this typically causes an 
additional cross-polarization. 

• Polarizer and OMT: a small part of the left circularly polar- 
ized signal could enter the right circularly polarized signal, 
and vice versa. This could heavily affect the V measurement 
while the effect on /, Q, U should be negligible. 

• LNA: a small fraction of the received radiation is reflected 
back to the previous components due to inevitable impedance 
mismatches, an effect called non-zero return loss. Once this 
spurious signal passes through the polarizer it faces the direc- 
tional coupler and is reflected back again. The two reflections 
added to a passage through the polarizer transform the circu- 
lar polarization into its orthogonal polarization. This leads to 
a spurious exchange of radiation between the left and right- 
hand channels. 

3.5. The down-conversion 

After the OMT, the two circularly polarized signals are ampli- 
fied and downconverted. The "conversion box" can be consid- 
ered linear (some deviation appears for very strong sources) but 
could introduce some polarization aberration: some small differ- 
ence in the signal path length could alter the phase difference 
between the two components, and the amplitude could also be 
sligh tly modified. The sign al, a fter the front-end according to 
IConwav & Kronbera 111 969(1 and lMcKinnonl Ill992ll . is then (see 
Table[T]for the adopted notation) 

_ g r (e r + d r e l ) 

V r = V ' (8a) 

_ g l (e l + d l e r ) 

V l = V ' (8b) 

V 1 + Z) i 

where G R = G R e~"^ R and Gl = Gl<? ^ l are the gain of the right 
and left channels and D R = D R e- itpR and D L = D L e~' l/ "- are the 
cross-talk factors. 

The receiver architecture usually offers narrow bandwidths 
(some tenths of MHz) at low frequencies (< 3 GHz) and wider 
bandwidths (some hundreds of MHz, few thousands in some new 
generation receivers) at higher frequency (>3 GHz). It can be 
shown that und er these circumstances the ba ndwidth depolariza- 
tion factor (see iGardner & Whiteoakl lll966ll ). which affects LP 
and is defined as the ratio between the observed degree of linear 
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polarization and the emitted one, is negligible for typical extra- 
galactic sources (a typical Rotation Measure of 50 rad/m 2 gives 
a depolarization factor of nearly 1 both at 2.5 GHz over a 100 
MHz band and at 5 GHz over a 500 MHz band). 

4. Correction of observations using the 
instrumental Muller matrix, T 

After the down-conversion, the signals are sent to the de- 
vices that perform the measurements by implementing Eqs. ©. 
Typical devices are square law detectors (to obtain the flux den- 
sities from the left and right channels) and multiplying polarime- 
ters (which supply Q and U, through a 90° shift of Vl or Vr). In 
App.|A]the calibration procedure is described (the complete cal- 
culations are detailed in Cenacchi (2009, PhD thesis, in prep.). 

4.1. Effects on the measurements 

The instrumental Muller matrix elements contain information 
about the amount of contamination among the four Stokes pa- 
rameters, due to spurious, unwanted, conversions arised inside 
the receiver. The conversion of Q, U and V into / is typically 
negligible, the conversion of / and V into Q and U represents 
a small but still significant corruption, the conversion of /, Q 
and U into V can completely overwhelm the result of the mea- 
surement, leading to the impossibility of disentangling the in- 
strumental contribution from the intrinsic V value. This clearly 
appears when dealing with completely unpolarized (thermal) 
sources. 

As an example, Table [3] lists the Stokes values measured in 
November 2007 during a test observation with the Effelsberg 
100-m telescope at 6 cm of two planetary nebulae. Being thermal 
sources they are expected to be unpolarized. The sources were 
observed with five beam wide on-the-fly cross-scans centred on 
the source. The Q and U values came directly from the multiply- 
ing polarimeter. The V value has been estimated as the difference 
between the left and right channels, after those were separately 
corrected for pointing (to compensate for the beam squint) and 
gain curve. The cross-scan technique intrinsically removes pos- 
sible differences between the internal channel noise (* 1 K for 
the Effelsberg telescope at 6 cm). 

The measured Q, U and V values, from unpolarized sources, 
are believed to be due to a contamination from Stokes /, that 
is, from a small net exchange in amplitude and phase between 
the left and right channels that, despite the multiplication and 
the correction for pointing and gain, still produces a spurious 
significant LP and CP detection. 

App.lBlcontains the definition of each element of the Muller 
matrix. 



4.2. How to apply the D-terms calibration 
4.2.1 . Derivation of the Muller matrix 

The proposed calibration scheme is completely based on ob- 
servation of astronomical calibration sources, each observation 
schedule must include a fully unpolarized source (Q = U — V — 
0) and a highly linearly polarized source (for which one does not 
require V = OK). The left and right channels must be corrected 
separately for pointing and gain curve. The raw V value can then 
be estimated as the difference of the measured kelvin. R and L 
opacity correction can be applied to the channels and the raw / 
can be derived as the sum of the R and L temperatures. The raw 
Q and U come directly from the multiplying polarimeter. 



By observing one or more unpolarized sources whose flux 
density is known, Eq. (|6]l is simplified as follows 



Im = tills 
Qm - hlh 
U m = hlh 
V m = ?4i/. s 



(9a) 
(9b) 
(9c) 
(9d) 



where tn,tzi,hi and ?4i are given in Eqs. ( IB. 11 1. These four 
matrix elements involve the parameters Di, Dr, ipi, <Pr, A c and 
*T, all of which must be known for the calibration to be applied. 

The polarization angle of the internal noise signal, A e is 
known (either from technical specifications or from previous 
measurements) and the remaining five quantities can be solved 
from six measurements on the calibrators (four Stokes parame- 
ters measurements on the unpolarized source and Stokes Q and 
U measurements on the linearly polarized calibrator). 

Once the D-terms, in amplitudes and phases, are determined, 
all terms of T can be calculated, as they depend on only these 
parameters. By applying this procedure, it is not necessary to 
make assumptions about the V value from the linearly polarized 
calibrator. 

4.2.2. Treatment of the errors 

Let the uncertainties measured on I m , Q„ u U m and V„, for the cal- 
ibrators) be a - /, ctq, o-jj and cry, then the corresponding uncer- 
tainty on the matrix elements are 



0"21 

o"3i 

0~41 



Is 



(TV 



h 

CTy 



(10a) 
(10b) 
(10c) 
(10d) 



Due to the non-linearity of the equations involved, it is not 
possible to propagate these errors through the estimate of the 
cross-terms, to the complete instrumental Muller matrix. It is 
anyway possible, by comparing the mathematical definitions of 
the matrix elements, to use the results ( TTOb to establish an upper 
limit for the remaining errors, e.g. it is clear from their definition 
that the error on mu and mi4 must be the same. By applying 
similar qualitative considerations, the complete error matrix can 
be estimated as 



( o"ll 3o- n cr n 

CT 2 i CT 2 i cr 2 i 0"2i 

cr 31 0"2i cr 2I cr 31 

cr 41 cr n 3o- n cr 41 ) 



(ID 



The intrinsic Stokes vector can be then obtained by inverting 
T, propagating step by step the errors cr nm during the inversion 
process and applying the following 



B^T ~ l S, 



(12) 



Finally, once S s is obtained, one or more flux-density cali- 
brators can be used to dermine the K/Jy conversion factor for I 
and V, and a strongly polarized calibrator can be used to deter- 
mine the K/Jy conversion factor for Q and U. 
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Table 3. Examples of raw measurements of all Stokes parameters for unpolarized sources at Effelsberg at 6 cm. Measurements were 
calibrated using a noise diode of known strength, but D-terms correction was not applied. 



Source 


L 


0"A 


Qm 


0"A 


u m 


0"A 


v m 


0"A 


LP 


o"a 


CP 


0"A 




[K] 




[K] 




[K] 




[K] 




[%] 




[%] 




NGC6572 


3.861 


0.008 


-0.016 


0.002 


0.006 


0.013 


-0.016 


0.004 


0.43 


0.12 


-0.42 


0.11 


NGC 7027 


17.100 


0.300 


-0.076 


0.011 


-0.014 


0.015 


-0.087 


0.022 


0.45 


0.07 


-0.51 


0.13 



5. Observed results 



5.3. Linear polarization 



5.1. Instrumental terms 

During the year 2007 several observations were carried out at the 
100-m Effelsberg telescope to test the new sub-reflector. Part of 
this test time has also been used to test the full Stokes polarimet- 
ric calibration. A selected sample of 43 sources (39 extragalac- 
tic sources, 2 planetary nebulae, 2 planets) have been observed 
nearly monthly at 6 cm. Subsequently, some observations at 1 1 
cm, 3.6 cm and 2.8 cm were also carried out. 

The Stokes parameters measurement listed in Table|4]for the 
planetary nebula NGC 7027 yielded also the following instru- 
mental terms (averages over time, and standard deviations) 

D L = (0.59 + 0.04)% 
D R = (0.16 + 0.08)% 
^ L = (-2.4 ±3.8)° 
cp R = (-3.7 ± 8.4)° 
y = (-4.1+2.7)°. 

These values are in good agreement with the receiver techni- 
cal specifications. The measured amplitudes are in good agree- 
ment also with EVN measurements. The definition of D-terms 
phase in the interferometry processing is quite obscure, so it is 
unclear whether the results are different or not. 



5.2. Circular polarization 

We used our D-term calibration to search for CP in 19 
AGNs with improved sensitivity over previous observations at 
Effelsberg. 

Table [5]lists a comparison at 6 cm between the raw circular 
polarization and the values obtained after the calibration, using 
NGC 7027 as unpolarized calibrator and 3C 286 as strongly po- 
larized calibrator. The data were observed in November 2007. 

After the D-term calibration, the second planetary neb- 
ula, NGC 6572, appears circularly unpolarized as expected. 
Significant amounts of CP were detected in the two known CP 
sources 0743-006 and 1519-273 at levels co nsistent with pub- 
lished levels (0743-006: V = (-0.51 +0.08)%.lAller et alj i200l 
and V = (-0.46 ± 0.05)%. iHoman et al.1 l200lll . both adopting 
the circular RL fr ame and for 1519-273: V = (-0.92 ± 0.17)% 
llAller et al.L[2003ll ). In the Effelsberg measurement of 15 19-273, 
the source appeared unpolarized from the raw values, but the 
D-terms calibration revealed a level of CP in good agreement 
with the values reported in the literature. The Effelsberg mea- 
surement of the weakly circularly polarized source 3C 279 also 
shows a level of C P close to that published ( V = (-0.17+0.01)% 
lAller et all 12001 ). 

T ime variabi l ity in CP and changes of sign are very com- 
mon llAlleretaUr2003ll . 1519-273 is the most stable CP source 
known so far. Table [6] shows a sample of CP values that we ob- 
served during 2007. Our detections also confirm the stability of 
the CP level of 1519-273. 



The proposed calibration scheme maintains the same reliabil- 
ity in the calibrat ion of linear polarization as that reported by 
iTurlo et all I1985H based on a 3x3 Miiller matrix. In addition 
our method supplies Stokes V, the D-terms and the complete er- 
ror treatment. 

The traditional polarization calibration is applied assuming 
a Miiller matrix unaffected by errors, whereas the scheme that 
we propose includes the error propagation also during the ma- 
trix inversion and its application, hence the resulting error esti- 
mates are intrinsically higher than those obtained with the stan- 
dard procedure. Tables [7] and [8] show the comparison between 
the standard and the new calibration schemes, when measuring 
linear polarization. 

The linear polarization values obtained are consistent with 
those obtained with the standard method. 



6. Conclusions 

A calibration procedure based on a 4x4 Miiller matrix has been 
proposed, that allows Stokes V measurements with single-dish 
telescopes that supply circularly polarized outputs. This method 
also provides Stokes /, Q and U, and so allows a full polarimetric 
study of sources under examination. 

The mathematical interpretation of the Miiller matrix ele- 
ments offers the possibility of deriving the instrumental polari- 
metric characteristics of the telescope (D-terms, in amplitude 
and phase) directly from the observation, allowing a direct com- 
parison with the technical specifications set by the engineers. 

The proposed procedure is valid for all single-dish telescopes 
equipped with a receiving architecture composed of a standard 
receiver, two total power detectors for the two circularly polar- 
ized channels and a multiplying polarimeter (one of the most 
common solution available nowadays), and can be applied at any 
observed frequency. 



Appendix A: The backend and the derivation of the 
measured Stokes parameters 

In the following, the average symbol <> is omitted, subscript s 
stands for source, and subscript c for calibration noise diode. 



A. 1 . Total power detectors 

To measure total power in the R and L polarizations, the mea- 
surements V R V* R = V 2 R and V L V* L = V 2 L are made. According 
to Eq. (O, introducing the square law detector gain (the scalar 
quantity g) and recalling that 



{XSf = A*B* 

(a + b)* = a* + zr 
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Table 4. Observed Stokes parameters for NGC7027 at 6 cm. Observations were made in 2007 at Effelsberg. Units are kelvin 



Month 


I 


cr 


Q 


cr 


u 


cr 


V 


cr 


February 


17.014 


0.143 


-0.079 


0.004 


-0.016 


0.006 


-0.106 


0.018 


April 


17.438 


0.289 


-0.070 


0.002 


-0.015 


0.005 


-0.133 


0.028 


May 


16.864 


0.079 


-0.073 


0.005 


-0.020 


0.008 


-0.108 


0.028 


June 


16.989 


0.061 


-0.066 


0.004 


-0.018 


0.008 


-0.122 


0.019 


August 


17.179 


0.245 


-0.066 


0.002 


-0.015 


0.007 


-0.124 


0.017 


November 


17.068 


0.343 


-0.076 


0.011 


-0.014 


0.015 


-0.087 


0.022 


December 


17.279 


0.265 


-0.075 


0.010 


-0.019 


0.015 


-0.093 


0.025 



Table 5. Observed circular polarization, before and after applying the D-terms calibration 



Source 




T 


cr 


1 / 

Vm 


(T 


J 

h 


cr 


1/ 

V S 


CT 






run 




|%J 








l%J 




aaca nn 




a 1 on 
4. loU 


A A 1 C 


A AA 


A 1 C 


A 1 AQ 

4. iyj 


A AO 1 


A TO 
U.ZO 


A 1 1 


A 1 1 /I i 

U134+3Z 


4o 


n no 
1 /.Zlo 


A 1 OQ 

U. loo 


A 1 A 


A 1 A 


1 H TOO 

1 /.Zoo 


U.zoo 


A 1 A 

U.Z4 


A AA 


A1 CO . -7 A 

Ul j 3+ /4 




J. 134 


a An 


A AQ 


A OA 

U.zU 


1 1 A H 

5. 14 / 


a 


A /II 

U.4Z 


A A/1 


\JJ lUTt 1 


R4 

. ' V... OH - 


df\ 4S6 


47R 

Ki.'-r 1 o 


-\).'-r'-r 


1 1 

U.l 1 


4fS f\S\ 


U.UJJ 


u.uu 


0f\ 


0212+73 




12.913 


0.120 


-0.36 


0.11 


12.967 


0.146 


0.08 


0.04 


0607-15 




11.352 


0.148 


-0.16 


0.07 


11.398 


0.201 


0.28 


0.08 


0716+71 




2.991 


0.028 


-0.28 


0.16 


3.004 


0.035 


0.18 


0.04 


0743-00 




6.063 


0.015 


-0.09 


0.10 


6.089 


0.132 


0.36 


0.08 


0835+58 




1.867 


0.027 


-0.21 


0.25 


1.880 


0.031 


0.25 


0.12 


0836+71 




6.692 


0.069 


-0.37 


0.15 


6.720 


0.085 


0.09 


0.08 


0951+69 


3C231 


10.320 


0.128 


-0.20 


0.15 


10.362 


0.117 


0.25 


0.04 


1226+02 


3C273 


117.944 


1.083 


-0.43 


0.09 


118.458 


2.049 


0.04 


0.09 


1253-05 


3C 279 


36.672 


0.223 


-0.33 


0.11 


36.831 


0.635 


0.11 


0.07 


1328+30 


3C286 


23.331 


0.305 


-0.28 


0.11 


23.428 


0.307 


0.14 


0.11 


1409+52 


3C295 


20.688 


0.439 


-0.20 


0.11 


20.763 


0.357 


0.25 


0.07 


1519-27 




4.700 


0.079 


0.20 


0.21 


4.719 


0.114 


0.62 


0.15 


1749+09 




10.873 


0.097 


-0.37 


0.12 


10.919 


0.238 


0.07 


0.10 


1800+44 




3.830 


0.079 


-0.49 


0.12 


3.844 


0.044 


-0.04 


0.05 


2134+00 




31.723 


0.672 


-0.45 


0.07 


31.836 


0.551 


-0.01 


0.07 


2145+06 




18.925 


0.211 


-0.39 


0.07 


19.003 


0.284 


0.05 


0.06 


NGC6572 




3.861 


0.008 


-0.42 


0.11 


3.878 


0.117 


0.03 


0.12 


NGC7027 




17.068 


0.343 


-0.51 


0.13 


17.130 


0.218 


-0.06 


0.05 


SATURN 




1.962 


0.028 


-0.85 


0.19 


2.165 


0.042 


-0.30 


0.09 



Table 6. Repeated observations of circular polarization, in percentage, from a sample of sources, during 2007. 





Feb 




Apr 




May 




Jun 




Aug 




Nov 




Dec 




Source 


V 


cr 


V 


cr 


V 


<x 


V 


cr 


V 


<T 


V 


cr 


V 


cr 


0743-00 


0.53 


0.04 


0.42 


0.05 


0.24 


0.08 


0.18 


0.06 


0.52 


0.04 


0.36 


0.08 


0.38 


0.06 


1253-05 


3C279 0.11 


0.05 


-0.12 


0.08 


-0.40 


0.09 


-0.25 


0.06 


-0.35 


0.10 


0.11 


0.07 


0.29 


0.12 


1519-27 






0.65 


0.14 


0.53 


0.08 


0.61 


0.11 






0.62 


0.15 






NGC6572 






-0.33 


0.19 


-0.06 


0.04 










0.03 


0.12 


0.19 


0.13 



one obtains 



*Rc 

V 2 I 

Vf 2 



The Teal is usually fully linearly polarized and character- 
ized by V c = OK. This simplifies the involved calculations. By 
comparing the measured signal with that of the noise diode, the 
receiver gain fluctuations are removed. In addition, the Teal cal- 
ibration procedure calibrates the source strength measurements 
in terms of antenna temperatures, in kelvin. Discarding the terms 
of order > 3 one obtains 



V R V* R = [El + D\E\ + E R E* L D* R + E L E' R D R ] (A. 1 a) 

1 + D R L J 

V L Vl = [E\ + D 2 L E 2 R + E L E R D* L + E R E* L D L ] (A. lb) 

1 + Di 1 1 



Eq. ( IA.U can be expressed in terms of the Stokes parameters 
or in terms of source linearly polarized flux Ii ps and phase dif- 
ferences. When the Teal calibration is applied, a common post- 
processing procedure allows one to transform the voltage at the 
output of the total power devices (TP) into kelvin, as follows 
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Table 7. Observed linear polarization, in percentage, from a sample of sources measured at 6 cm at Effelsberg during 2007. 
Superscripts 1 indicates the results obtained with the new calibration procedure, 2 those obtained with the traditional calibration 
procedure. 







Feb 




Apr 




May 




Jun 




Aug 




Nov 




Dec 




Source 




P 


cr 


P 


cr 


P 


er 


P 


cr 


P 


O" 


P 


cr 


P 


O" 


0743-00 1 




0.81 


0.19 


0.87 


0.12 


0.97 


0.54 


1.28 


0.27 


1.12 


0.17 


1.11 


0.60 


0.84 


0.53 


0743-00 2 




0.82 


0.13 


0.85 


0.09 


0.97 


0.20 


1.24 


0.06 


1.06 


0.08 


0.98 


0.09 


0.90 


0.13 


1253-05 1 


3C279 


3.25 


0.28 


1.24 


0.21 


0.81 


0.45 


0.96 


0.40 


1.14 


0.36 


1.42 


0.86 


1.11 


1.63 


1253-05 2 


3C279 


3.15 


0.08 


1.17 


0.03 


0.82 


0.03 


0.91 


0.02 


1.11 


0.04 


1.16 


0.04 


1.06 


0.00 


1519-27 1 








3.62 


0.71 


2.96 


0.25 


5.79 


0.58 






4.70 


0.62 






1519-27 2 








3.50 


0.01 


3.20 


0.16 


5.83 


0.01 






4.55 


0.27 






NGC 6572 1 








0.09 


0.79 














0.21 


1.27 


0.31 


0.77 


NGC 6572 2 








0.16 


0.01 














0.24 


0.27 


0.48 


0.59 



Table 8. Observed polarization angle, in degrees, from a sample of sources measured at 6 cm at Effelsberg during 2007, superscripts 
1 indicate the results obtained with the new calibration procedure, 2 those obtained with the traditional calibration procedure 







Feb 




Apr 




May 




Jun 




Aug 




Nov 




Dec 




Source 




PA 


O" 


PA 


cr 


PA 


cr 


PA 


cr 


PA 


cr 


PA 


cr 


PA 


cr 


0743-00 1 




29.5 


0.4 


31.4 


0.2 


32.2 


0.3 


28.1 


0.4 


30.6 


0.2 


31.0 


1.0 


26.0 


0.6 


0743-00 2 




28.8 


1.8 


30.6 


2.1 


31.8 


1.4 


29.6 


1.4 


29.8 


1.2 


34.4 


2.3 


24.6 


7.5 


1253-05 1 


3C 279 


-28.9 


0.1 


-28.7 


0.2 


-36.5 


0.4 


-17.7 


0.4 


-1.3 


0.4 


12.0 


0.5 


9.6 


0.8 


1253-05 2 


3C 279 


-29.9 


0.6 


-28.8 


1.3 


-31.6 


1.7 


-17.9 


0.6 


-1.9 


1.0 


13.5 


0.5 


8.6 


0.1 


1519-27 1 








3.3 


0.1 


-34.3 


0.1 


-36.4 


0.2 






130.1 


0.3 






1519-27 2 








2.2 


0.3 


-28.1 


6.6 


-36.6 


0.1 






129.4 


1.8 







7 '« = 7T77 — \l s [l+D 2 R + D 2 L + D L cos(cc L + A c ) 
UN iv 

+D R cos((p R - Ac)] 

+Qs [Dr cos ip R + D L cos ip L + 2D R D L cos tp R cos{tp L + A c ) 
+2D R D L cos <p L cos((p R - Ac)] 

+ U s [D R sin ip R + 2D R D L sin <p R cos(<^ L + A c ) - D L sin <p L 
-2D R D L sin <p L cos(ip R - A c )] 

+V S [D 2 R -D\-D L cos((p L + Ac) + D R cos((f R - A c )]} 

(A.3a) 

Vm = 7T7r~ i J s[ D R cos(if R - A c ) - D L cos(cp L + A c )] 
UN iv 

+Qs [Dl cos ip L + 2D R D L cos co L cos(co s - A c ) - D R cos <p R 
-2D R D L cos cp R cos((f L + A c )] 

+ U S [-D L sin <p L - 2D R D L sin cp L cos(^ - A c ) - D R sin ip R 
-2D R D L sin (p R cos(cp l + A c )] 

+V s [l + D R cos((p R - Ac) + D L cos(cp L + A c )]} (A. 3b) 

where 

DN IV = 1 + D R + u\ + 2D R cos((p R - A c ) + 2D L cos(<p L + A c ) 
+ AD R D L cos((p R - A c ) cos(cc L + A f ). 

These two equations show how the measured /,„ and V m are 
contaminated by I s , Q s , U s and V s depending on the values of 
the D-terms. 



A.2. Multiplying polarimeter 

To measure the linear polarization components, Q m and U,„, the 
measurements V R V* L + VlV r and i (y R V* L - Vl,V r ^J are made using 
a multiplying polarimeter. The outputs of the multipliers in the 
polarimeter are related to the incoming electric fields by 



VrVI = [ ErE * l + DrE2l + D ' lE2r + DrD l ElE k\ (A ' 4a) 

V l V *R = 7^ [ElE'r + Dr E 1 + D L E 2 R + D' R D L E* L E R ] (A.4b) 



DN M p 
where 

DN M p = sfl + D\ + D\ + D 2 R D 



The Eqs. JA.4| > can be expressed in terms of Stokes parameters or 
in terms of linearly polarized flux density. By assuming that gQ = 
gu and jq - yu (standard procedure at Effelsberg), it is possible 
to apply to both the Q and U channels the same calibration factor, 
according to the following 



(A.5) 



where Q c and U c are the Teal signals coming through the Q and 
U channels and measured at the end of the backend. 



The Teal is applied analogously to Eqs. (IA.2t . The ratio be- 
tween the Ii pc value and the corresponding measured voltage 
gives the conversion factor K/V to be applied to the on-source Q 
and U measurements 



Qm = 
U m = 



<Q C 2 + U 2 

U s 



Hpc 



• llpc 



(A.6a) 
(A.6b) 



Ve, 2 + u 2 

Discarding the terms of order > 3 the following is obtained 
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• Propagation of Q s into Q m 



Q„ 



1 



DN, 



QU 



■ [I s [D L cos (*P + <p L ) + D R cos OF - <p R )] 



+Q S [cos ¥ + D R D L cos OF + <£)] 

+U S [sin ¥ - D R D L sin OF + O)] 

+V S [D s cos OF - <Pr) - D L cos (¥ + <p L )]} 



(A.7a) 



U m =- 



1 



■ {/, [-D L sin OF + <Pl) ~ D R sin OF - <f R )] 



+Q S (- sinT - DrDl sin OF + O)] 

+ t/ 5 (cos ¥ - D R D L cos OF + <D)] 

+y. s (-Dr sin OF - + D L sin OF + <p L )]] 



(A.7b) 



where 



DNqu = [l + £>l + £>r + 2Z) L £> S cos (<D + A) cos A 

1 

+2D, cos (ifii + A) + 2Dr cos (A — (p R )] 2 . 

The channels Q and U could also be calibrated separately. In 
this case in Eqs. ( IA.7t two different denominators, one for each 
channel, would be present. 

With Eqs. dA.3b and ( IA.7I ) we have all 16 elements of the 
Miiller matrix in terms of D-terms, required to relate measured 
and true Stokes parameters. 

Appendix B: Definitions of the matrix elements 

From Eqs. dA.3t . ( lA.7b and recalling the definition (0, the coef- 
ficients of the instrumental Miiller matrix T can be summarized 
as follows: 

• Propagation of I s into /,„ 

f n = 7TT7 — [ 1 + D \ + D l + °l cosCp L + A c ) + D R cos(ip R - A c )l 
UN iv 

(B.la) 



Propagation of Q s into /„, 
1 



f i2 = [Dr cos 9? s + Di cos yp L + 2D R D L cos ^ cos(y> L + A c ) 

(B.lb) 



DNiv 

+2D R D L cos tp L cos((fi R - A c )] 



Propagation of U s into I m 
1 



f i3 =t^— [A? sin ^R + 2D R D L sin ^ s cos(^ L + A c ) - D L sin <p L 

(B.lc) 



-2D R D L sin </> L cos(<£> s - A e )] 



Propagation of V s into /, 
1 



fl4 



[D s -D 2 l -D l cos(ip L + A c ) + D R cosiifR - Ac)] 

(B.ld) 



Propagation of / s into 2, 
1 



*M = 



DN, 



QU 



[D L cos OF + <p L ) + D R cos OF - <p R )] 



(B.le) 



ti_2 



1 



[cos T + D R D L cos OF + <D)] 



DN QU 

Propagation of U s into Q 



?23 = 



1 



[sin^-DfiDi sin OF + O)] 



DN QU 

Propagation of V s into Q, 



tu - 



1 



[D R cos OF - <p R ) - D L cos OF + <p L )] 



DNqu 
Propagation of I s into U, 
1 



*31 = 



DN, 



QU 



[—Dl sin OF + <p L ) - D R sin OF - tp R )] 



< Propagation of Q s into U„ 
1 



hi 



DN 



QU 



[- sin "¥-D R D L sin OF + <D)] 



hi 



Propagation of U s into U„ 
1 



[cos T - D R D L cos OF + <£)] 



ZW e£ / 

Propagation of V. s into I/, 
1 



f 34 - 



[-D R sin OF - ^) + D L sin OF + <f L )] 



DNqu 
Propagation of I s into V, 
1 



?41 - 



[D R cos(<p R - A c ) - D L cos(<p L + A c )] 



Propagation of Q s into V, 
1 



(B.lf) 



(B.lg) 



(B.lh) 



(B.li) 



(B.lj) 



(B.Ik) 



(B.ll) 



(B.lm) 



f 42 =7^— [Dl cos ^ l + 2D R D L cos <p £ cos(^ - A e ) - D R cos ^ 

(B.ln) 



DN IV 

-2D R D L cos ^p fi cos(^ L + A c )] 



Propagation of U s into V, 
1 



f 43 = -^-rr~ [~ d l sin ^ £ - 2D R D L sin 9? L cos(<#? - A C )-D R sin 

(B.lo) 



DA^/y 

-2D R D L sin </3 S cos(<^l + A c )] 



Propagation of V s into V, 
1 



f44 - 



DAT/y 



[1 + D R cos(tp R - A c ) + D L cos(<p L + A c )] (B.lp) 
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